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The light vector mesons (p, u, and <f>) were produced in deuterium, carbon, titanium, and iron 
targets in a search for possible in-medium modifications to the properties of the p meson at normal 
nuclear densities and zero temperature. The vector mesons were detected with the CEBAF Large 
Acceptance Spectrometer (CLAS) via their decays to e + e~ . The rare leptonic decay was chosen 
to reduce final-state interactions. A combinatorial background was subtracted from the invariant 
mass spectra using a well-established event-mixing technique. The p meson mass spectrum was 
extracted after the w and (f> signals were removed in a nearly model-independent way. Comparisons 
were made between the p mass spectra from the heavy targets (A > 2) with the mass spectrum 
extracted from the deuterium target. With respect to the p-meson mass, we obtain a small shift 
compatible with zero. Also, we measure widths consistent with standard nuclear many-body effects 
such as collisional broadening and Fermi motion. 
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I. INTRODUCTION 

Quantum chromodynamics (QCD), the theory of the 
strong interaction, has been remarkably successful in 
describing high-energy and short-distance-scale experi- 
ments involving quarks and gluons. A major difficulty 
for QCD has been its application to low-energy and large- 
distance-scale experiments. However, the symmetries of 
QCD (such as chiral symmetry) provide guiding princi- 
ples to deal with strong interaction phenomena in the 
non-perturbative regime, where the strength of the in- 
teraction increases quickly. 

Much of the hadron mass is generated dynamically, 
and these masses are somewhat effected by the sponta- 
neous breaking of chiral symmetry. For example, the 
proton has a mass of approximately 1 GeV, that is 
much larger than the summed masses of its constituent 
quarks, which are a few MeV. In the early 1990's, various 
models 0, 0] related the in-medium hadron properties 
with chiral symmetry restoration. At high temperatures 
and/or densities, chiral symmetry is restored as the chiral 
condensate, (qq) , approaches zero Q . A partial restora- 
tion of the chiral symmetry is predicted to be achieved 
at normal nuclear densities, such as inside a heavy nu- 
cleus Q. 

In this article, we study the in-medium properties 



(mass and width) of the p meson. Our experiment is 
the photoproduction of the light vector mesons in nu- 
clei and their decay into e + e~ pairs with a search for 
medium modifications of the p meson. This experiment 
is the first to use an electromagnetic interaction in both 
the production and decay channels. The incident beam 
was a tagged photon beam, which has the advantage of 
interacting through the entire nuclear volume. Hadronic 
beams, such as protons and pions, interact on the surface 
of the nucleus. The e + e~ pairs are preferable decay chan- 
nel over hadronic channels. Hadronic final-state interac- 
tions will distort the reconstructed mass spectrum and 
mask any signals from chiral symmetry restoration. The 
leptonic decay eliminates final-state interactions. This 
experiment was conducted with the CEBAF Large Ac- 
ceptance Spectrometer (CLAS) at Jefferson Lab. The 
CLAS detector, with its excellent electron and positron 
identification, is ideal for this measurement. Our results 
have recently been published in Physical Review Let- 
ters @ . This article will describe in greater detail the ex- 
perimental techniques and analysis procedures employed 
to arrive at our result. 

The first experimental results on p-meson medium 
modifications came from relativistic heavy-ion experi- 
ments. The CERES @ collaboration at CERN reported 
an excess in the e + e™ mass spectrum in the p-meson 
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region. The experiment was a comparison of the mass 
spectrum from the p-Au reaction to the mass spectrum 
from the Pb-Au reaction. The data from the proton- 
induced reaction were described well by incorporating 
a cocktail of hadronic decay channels into their analy- 
sis. In the CERES work, the heavy-ion collision data 
displayed an enhancement in the mass range between 
300 and 700 MeV. This result could be explained as a 
temperature-induced decrease in the mass of the p me- 
son [3] • A second CERES measurement [|| with improved 
mass resolution confirmed the previous result. Also at 
the CERN-SPS, the HELIOS/3 collaboration studied 
the di-muon mass spectrum up to the J/\l/ mass with 
proton and sulfur beams on a tungsten target. They ob- 
served an excess in the di-muon mass spectrum below the 
</>-meson mass with the S-W reaction as compared with 
the p-W reaction. Moreover, the NA60 collaboration 
reported a doubling of the p-meson width from their di- 
muon measurement from In-In collisions with no change 
in the p mass. This result appeared to confirm the pre- 
diction of Ref . [ill. Recent publications from the NA60 
collaboration |l2l. Il3|| confirm the earlier result [icj |. 

Interpretation of relativistic heavy-ion collisions is 
complicated because the reaction occurs in a non- 
equilibrium state before proceeding to equilibrium, while 
theoretical models predict the hadronic properties at 
equilibrium (normal nuclear density and zero tempera- 
ture). The medium modifications of the p, w, and <p 
mesons are predicted to be large enough to be observed 
in elementary reactions with hadron and photon beams. 
Bertin and Guichon [l4[ predict a 120 MeV shift in the 
p-meson mass from their quark model of the nucleus. 
Brown and Rho [l[ constructed an effective chiral La- 
grangian with suitable QCD scaling and calculated a 20% 
decrease in the mass of the p meson. Hatsuda and Lee [2| 
used QCD sum rules to obtain the following parameteri- 
zation for the masses of the three light vector mesons: 

m = mo ( 1 — a— ) , (1) 
V Pa J 

where m is the mass in the medium, mo is the mass in 
the vacuum, p is the nuclear density, and po is the nor- 
mal nuclear density. The variable a parameterizes the 
mass shift. The current prediction of a from Ref. [2[ 
is 0.16 ± 0.06. In the past decade, more sophisticated 
models have been developed. However, some calcula- 
tions predict a decrease in the mass fla. Hll Ht} . oth- 
ers claim an increase [II [H, EE ED, EaEJ and some 
even predict the appearance of multiple peaks in the mass 
spectrum [24l EH- In addition, nuclear many-body ef- 
fects have been incorporated in these models. Ref. [llj 
predicts a doubling of the width of the p meson in a 
hot, dense environment, such as produced in relativistic 
heavy-ion collisions. 

The TAGX collaboration reported a large decrease of 
the p-meson mass in the reaction 3 He(7, 7r + 7r~)X, where 
the pion pairs result from sub-threshold p-meson produc- 
tion and decay 26]. Their tt + tt— mass spectrum was fit 



with "non-p" processes, a medium-modified p process, 
and a p process that is not modified by the medium. 
The best fits were consistent with the 20% reduction pre- 
dicted by Brown-Rho scaling. A second analysis based 
on the longitudinal polarization of the p mesons reported 
a smaller but significant decrease in mass [27| . These 
results are questionable given the small density of the 
nucleus and the final-state interactions on the pion pairs. 

The KEK-PS collaboration detected p, uj, and <f> 
mesons decaying into e + e _ pairs from a reaction of 12- 
GeV protons incident on C and Cu targets [H, ES E3] • 
They reported a decrease in the p-meson mass consistent 
with a = 0.092±0.002. The Crystal Barrel/TAPS collab- 
oration reported a decrease in the mass of the u) meson 
by studying the photoproduction of low-momentum u> 
mesons in Nb that decayed through the tt°j channel [3l| . 
The decrease in mass translates into a on the order of 
0.13 for an average density of 0.6po- 

The article is organized in the following manner. Sec- 
tions [TT] and IIIII describe the experimental set-up and 
the lepton pair identification, respectively. In Section [IV! 
background contributions are discussed. Section [V] re- 
views the realistic generator and simulations employed in 
extracting the p mass spectrum. SectionslVTlandlVTIlpro- 
vide the final results and a study of the systematic errors. 
In Section IVIII1 our results are compared quantitatively 
with theoretical calculations and other experiments. 



II. EXPERIMENTAL DETAILS 
A. Detector Components 

This experiment was conducted in Hall B at the 
Thomas Jefferson National Accelerator Facility (TJ- 
NAF). An electron beam that was accelerated by the 
Continuous Electron Beam Accelerator Facility (CE- 
BAF) [32] to an energy of 3.062 or 4.023 GeV was di- 
rected into Hall B. The RF timing structure of the elec- 
tron beam is in 2 ns bunches. The beam was incident on 
a radiator with a 10 -4 radiation length thickness and de- 
flected into the Hall B photon tagging facility [33]. The 
outgoing bremsstrahlung photon beam was reduced with 
a 1 mm-aperture collimator such that the beam-spot size 
on the target was 1 cm. The incident tagged-photon 
flux on target was approximately 5 x 10 7 photons/s over 
an energy range from 0.61 GeV to 3.82 GeV. A multi- 
segmented target was used and will be discussed in the 
following section. The outgoing vector mesons were re- 
constructed from their decay to e + e~ pairs. The lep- 
tons were detected with the Electromagnetic Calorime- 
ters (EC) [13 and Cerenkov Counters (CC) [35] of the 
CLAS detector [H in Hall B. 

The CLAS detector is ideal for this experiment owing 
to its ability to detect multi-particle final states and its 
high e/ir rejection factor (see Sec. IIII B[) . The primary 
components of the spectrometer are a six-coil supercon- 
ducting toroidal magnet, a series of 3 Drift Chamber 
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modules (DC) [37|, time-of-flight Scintillation Counters 
(SC) [Hi, the EC, and the CC. The magnet produces a 
mostly azimuthal magnetic field around the middle re- 
gion of the DC. The momentum resolution is 0.5-1% for 
charged particles, depending on the kinematics. To re- 
duce the low-energy e~ and e + background from pair 
production in the targets, a "mini-torus" magnet was sit- 
uated just beyond the target region and inside of the DC. 
The CLAS detector is divided into six identical spectrom- 
eters that are oriented in the 9 direction and are refered 
to as sectors. 



B. The Target 

The purpose of this measurement was to produce vec- 
tor mesons from a variety of target nuclei. The tar- 
get materials were liquid deuterium (LD 2 ), carbon, ti- 
tanium, iron, and lead. To reduce systematic errors due 
to beam intensity and position, the target assembly was 
constructed such that each material was simultaneously 
in the beam. Moreover, the carbon target was divided 
into four foils and interspersed between the heavier tar- 
gets to study systematic effects of target position. The 
order of target materials from upstream to downstream 
was LD 2 , C, Fe, C, Pb, C, Ti, and C . The diameter of 
each target foil was 1.2 cm to match the 1 cm beam-spot 
size. The distance between each target foil was 2.5 cm. 
The LD 2 cell was 6.2 cm in length with a 1.64 cm diame- 
ter entrance window. This spacing was optimized to min- 
imize any multiple scattering of the electron or positron 
in subsequent target foils. The total thicknesses of the 
LD2, C, and Pb targets were approximately 1 g/cm 2 , 
and the Ti and Fe foils were about 0.5 g/cm 2 each. In 
order to increase the statistics in the present analysis, 
the Ti and Fe data were combined. This is justified since 
their radii (w A 1 / 3 ) are similar. Each target type had 
about the same number of nucleons. Due to their narrow 
widths, the to and <fr mesons were used to normalize the- 
oretical mass-shapes and were subtracted from the total 
mass spectrum. 

The reconstructed vertex position of the e + e~ pair was 
employed to determine in which target the vector meson 
was produced. Fig. [1] shows the reconstructed vertex po- 
sition from the data. The vertex resolution of CLAS is 
enough to clearly identify each target. 



III. EVENT SELECTION 

A number of cuts were applied to obtain the clean- 
est possible sample of lepton pairs and to remove any 
hadronic backgrounds. The cuts were grouped into two 
categories. The first set of cuts employed the EC and 
the CC for lepton identification and was applied to each 
particle individually. The second set of cuts matched 
the particle pairs by their interaction vertex position and 
time. Before any cuts were applied, the data were filtered 




z (cm) 

FIG. 1: (color online) Reconstructed e + e~ vertex position 
along the beam direction, z. The position of z — corre- 
sponds to the center of the CLAS detector. 



to select events with two oppositely-charged particles in 
CLAS. 



A. Lepton Identification 

Since almost 100% of the p mesons decay into two pi- 
ons and the relative branching ratio for the e + e~ decay 
channel is of the order of 10 -5 , it was crucial to discrim- 
inate between e + e~ and 7r + 7r~ pairs. In the EC, the lep- 
tons produce electromagnetic showers and deposit energy 
that is proportional to their momentum (see Figs. and 
[3]) with small corrections for momenta below 1 GeV. The 
pions are minimum-ionizing in the EC and deposit a con- 
stant amount of energy per unit length. A momentum- 
dependent cut was applied in EC deposited energy to se- 
lect the lepton band. The EC is segmented longitudinally 
into two sections. The EC inner section has fewer layers 
than the outer section; therefore, the minimum-ionizing 
pions lost the majority of its energy in the outer section. 
A cut was applied to exclude particles that deposited less 
than 45 MeV in the inner EC section. The lepton iden- 
tification was refined further with a coincidence between 
the EC and CC. 

To test the level of pion contamination, the EC/CC 
cuts were applied to a well-identified sample of pions. 
The sample chosen was ir~ from the well-known A decay 
in the reaction jp — > KK + — ► pir~ K + . First, a por- 
tion of the data was filtered for a final state containing 
a proton, a 7r _ , and a K + . The second step was to re- 
construct the pir~ invariant mass (see Fig SI) • The last 
step was to apply the EC and CC cuts to the pions in 
the mass range from 1.11 to 1.12 GeV (i.e. all the pions 
that contribute to the A peak). A sample of 5581 "good 
pion" candidates was selected in this fashion. From the 
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FIG. 2: (color online) The deposited energy in the EC versus 
momentum for positively-charged particles. 
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FIG. 4: (color online) Invariant mass of pn~ from the reaction 
7P — * p-K~ K + . The shaded region is the cut to select the A 
events. 




FIG. 3: (color online) The deposited energy in the EC versus 
momentum for negatively-charged particles. 



invariant mass plot, there was a small amount of back- 
ground under the peak. Of those 5581 7r~ particles, 3 
passed the EC and CC electron identification cuts. The 
upper limit of the pion rejection factor was estimated to 
be 5.4 x 1CP 4 . The analysis requires a lepton pair in the 
final state, and this translates into a rejection factor of 
2.9 x 10~ 7 for the pair. The pion momenta from this 
AK + test were less than 0.5 GeV where the lepton and 
pion bands have the maximum overlap. Thus, this anal- 
ysis sets an impressive upper limit in the most difficult 
kinematic range for lepton identification. 

An additional test was done to study the effect of the 
lepton identification cuts on n + and 7r _ particles with 
high momentum. Single charged-pion events were ana- 
lyzed, and the number of tt + and tt~ (identified by par- 



ticle identification from the SC) were counted before and 
after applying the EC/CC cuts. The ratio of the number 
of surviving 7t + and 7r~ (misidentified as leptons) after 
the cuts to the number of ir + and tt~ before the cuts was 
9.5 x 10~ 5 and 5.7 x 1CT 4 , respectively. 



B. e e Pair Identification 

In order to identify the lepton pairs from the decay of 
vector mesons, a set of cuts was applied on the interac- 
tion vertex. The e~ and e + were matched in position 
along the beam direction (z), radial position in the tar- 
get, and production times. The production time of each 
particle is measured from its arrival time at the SC and 
its trajectory to the vertex. 

Because the lepton pairs from the decay of the vec- 
tor mesons have large opening angles, we reduced the 
background from competing processes by requiring that 
the two leptons be detected in different sectors of the 
CLAS detector. This requirement removed the large 
background from the pairs produced by pair production, 
Bethe-Heitler processes, and tt° and rj Dalitz decays that 
have a small opening angle. The different sector cut re- 
moved about 40% of the events that survived the vertex 
and timing cuts. The e + e~ invariant mass is shown in 
the top panel of Fig. [5] for all targets and before apply- 
ing any cuts on the lepton pair candidates. The middle 
panel of the same figure shows the invariant mass dis- 
tribution after the vertex position and timing cuts, over- 
laid with the distribution of the same-sector e + e~ events. 
One final cut was to exclude events with low-momentum 
leptons. The low-mass background is dominated by lep- 
tons with momenta below 0.5 GeV. The bottom panel of 
Fig.[5]shows the e + e~ invariant mass after all of the cuts. 
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The background with a mass below 0.5 GeV was greatly 
reduced after the low- momentum cut, while the higher 
mass region was not affected. The low-momentum cut 
also produced similar e~ and e + acceptances. The mo- 
mentum of the leptons was corrected to account for minor 
misalignments in the toroidal magnet and for energy loss 
in the target materials. It was difficult to extract the p- 
meson signal from Pb due to high backgrounds and the 
absorption of the u> and 4> mesons, so the Pb data were 
discarded. 
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IV. BACKGROUND STUDIES 

Possible background channels in the invariant mass 
spectrum were investigated. The background may be a 
mixture of same photon (correlated) and different photon 
(uncorrelated) interactions that produce e + and e~ . The 
probability of an untagged photon and a tagged photon 
being in the same RF timing bunch is about 25%. Ex- 
amples of correlated sources of e + e~ pairs, in addition 
to p, w, and <f> mesons, would be uj and 77 Dalitz decays 
[u> — ► 7r°e + e~ and 77 — ► e + e~j). These channels were 
simulated with our physics model described in Sec. [V] 
There we discuss the possibility of contributions from 
Bethe-Heitler processes and 7r° Dalitz decay to the e + e~ 
mass spectrum. These channels were found to have negli- 
gible contributions to the invariant mass spectrum above 
0.5 GeV. Uncorrelated background is caused by incident 
photons from out-of-time beam bunches. An example 
would be an electron from a photon conversion com- 
bined with a positron from a vector-meson decay. The 
model-independent determination of the shape and nor- 
malization of the uncorrelated background will also be 
discussed. 
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FIG. 5: (color online) The e + e~ invariant mass spectrum 
summed over all targets before applying any cuts to pair can- 
didates (top), after vertex position and timing cuts (middle), 
and after the vertex position, timing, lepton momenta, and 
different-sector cuts (bottom). The overlaid histogram in the 
middle plot shows events with the e~ and e + detected in the 
same sector of CLAS. 



A. Bethe-Heitler Processes 

Elementary Bethe-Heitler processes are shown in 
Fig. [H The incoming photon produces a lepton pair, and 
one of the two leptons interacts with a target nucleon via 
exchange of a virtual photon. Although the majority of 
these pairs have very small opening angles and were cut 
out from our data sample by requiring that the electron 
and positron were detected in different CLAS sectors (as 
discussed in Sec. MI Bp . a full study was completed to 
account for any remaining events after the cut. A Bethe- 
Heitler event generator was developed that incorporated 
the experimental conditions with the CLAS detector. 





FIG 6: Feynman diagrams for the Bethe-Heitler processes. 
An incoming photon with four-momentum k produces a lep- 
ton pair. The negative and positive leptons have four momen- 
tum p_ and p+ . The initial and final momenta of the nucleon 
are denoted Pi (p/). The four-momentum which is exchanged 
between the nucleon and one of the leptons is denoted by q. 

We recently calculated the Bethe-Heitler cross section 
on a nucleus. The momentum distribution of the nu- 
cleons inside the nucleus was described by a Fermi gas 
model, and a Woods-Saxon distribution was employed to 
approximate the nuclear density distribution. 
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In order to translate the cross section into the expected 
number of events, the CLAS-accepted events were nor- 
malized to the number of target nuclei and the number 
of incident photons as a function of beam energy and 
weighted by the Bethe-Heitler cross section. The simu- 
lation showed that Bethe-Heitler processes made a neg- 
ligible contribution to the experimental yield (less than 
0.01% in the region of the p meson in the invariant mass 
spectrum for a Fe target). 



B. Dalitz Decay of 27r°'s 

Fig. [7] shows two 7r°'s decaying via a Dalitz decay 
(7T° — > e + e~7) forming an e + e~ background. The Dalitz 




uncorrelated sources. The salient feature is that they pro- 
duce the same-charge lepton pairs as well as oppositely 
charged pairs. The same-charge pairs (e + e + and e~e~) 
provide a natural normalization of the uncorrelated back- 
ground. The CLAS detector measured lepton pairs of any 
charge combination, so the data contained both same- 
charge and opposite charge pairs. This combinatorial 
method has also been used in the past for measurements 
of opposite-sign pions and muons [4l], |42j • This method 
has also been used in the extraction of resonance sig- 
nals [43| and proton femtoscopy of eA interactions [44[ ■ 
The combinatorial background is statistically approx- 
imated by an event-mixing technique. An electron of a 
given event is combined with a positron of another event, 
as the two particles are completely uncorrelated. This 
procedure produces a sample of random e + e~ pairs that 
have the same trigger conditions and acceptance as the 
data. 

In this study we followed the prescription used by e + e _ 
measurements in Ref. [45[ . The assumptions in this tech- 
nique are that the two leptons have similar acceptances 
and that there are only two e + e~ pairs from uncorrelated 
sources. In general, the probability of having r leptons 
in an event with an expectation value of p, is given by 
a Poisson distribution, so that the single lepton event 
probability for the positron (P+) and electron (P-) are 
given by: 

P( r , M ) = ^l e -M ( 2 ) 



FIG. 7: Feynman diagrams of the Dalitz decays of two r's. 
The combination of an electron and a positron from separate 
decays (dashed oval) could lead to an e + e~ background. 

branching ratio is about 1% [HI, and the tt° tt° cross 
section is about 10% of the p cross section [4(|. Thus, 
the semi-correlated, e + e~ yield from the Dalitz decay 
of 27r°'s will be of the same order of magnitude as the 
p — > e + e~ channel. 

The 7T° 7T° events (jp — > 7r°A + — > 7r°7r°j>) were simu- 
lated by an exponential t distribution between the beam 
photon and one of the 7r°'s. The total cross section was 
taken from Ref. [40j and integrated over the total photon 
flux. The leptons were propagated through the CLAS 
simulation software, and 0.02% of the generated events 
survived the selection cuts. In the invariant mass range 
above 0.2 GeV, this background channel was found to be 
negligible. 



C. Combinatorial Background 

Experimentally, the true lepton pair production has a 
background of random combinations of pairs due to the 



P + =P(l,p + ) = p+e-^ (3) 

P-=P(l 1 p-) = P-e-^-. (4) 

The probability of two same charge leptons in an event 
is: 

P ++ =P(2,p + ) = ^ 2 + e-^ (5) 

P^=P{2,p_)= l -ple-»-, (6) 

and the probability of having two opposite-charge leptons 
in an event will be: 

P+- = P+P- (7) 

Replacing P + and P_ and rearranging gives for e + e + and 
e~e~ events: 

P++ = 5A*+e"" + = \p+e'^p +e ~^e^ = \p\e^ 

(8) 

P__ = ipV-, (9) 
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where, ps 10 3 <C 1 => e' 1 — » 1. Therefore, 

P+_ = 2y/ P++P— (10) 

Since these probabilities are proportional to the number 
of events, Eq. [TO] was used to obtain the actual number 
of opposite-sign pairs that contribute to the mixed-event 
background. 

The data were also filtered to require two and only two 
same-charge leptons in each event. The same-charge lep- 
tons in each event are then combined to obtain the e + e + 
and e~ e~ invariant mass spectra. These distributions 
give the normalization for the e + e~ combinatorial back- 
ground. In this analysis, the statistically smoothed shape 
of the combinatorial background was obtained by mix- 
ing single-lepton events many times. The single-lepton 
events are chosen as uncorrelated sources of opposite- 
charge leptons. The normalization of the combinatorial 
background for each target was then obtained by count- 
ing the number of same-charge pairs and calculating the 
number of expected opposite-charge pairs from Eq. 1101 
The single-lepton mixed-event distribution was then nor- 
malized to the number of expected opposite-charge pairs. 
The uncertainty of the normalization was estimated at 
±7%, limited by the number of the same-charge lepton 
pairs. 

The result is shown in Fig. [5] for the 2 H, C, and com- 
bined Fe — Ti targets. The combinatorial background 
describes the shape of the background very well and is 
self-normalizing. This is an excellent method compared 
to other analyses where normalization of the mixed-event 
background was obtained by fitting the mixed-event dis- 
tribution to the data [HHHlU. 



V. EXTRACTION OF THE p MASS SPECTRA 
A. Simulations 

To simulate the expected physics processes, a realis- 
tic model was employed and corrected for the CLAS ac- 
ceptance. The events were generated using a code based 
on a semi-classical Boltzmann-Uehling-Uhlenbeck (BUU) 
transport model [H, |4?| that was developed at the Uni- 
versity of Giessen. This model treats nuclear effects such 
as the shadowing of the photon-induced reactions, as 
well as modeling of the Fermi motion, Pauli blocking, 
Coulomb interaction, final state interactions, and colli- 
sional broadening. Additional medium effects, such as 
decreasing the mass of the vector mesons [2|, can be in- 
corporated on demand. In photoproduction reactions, 
the meson is formed throughout the nuclear volume, un- 
like reactions with hadronic beams, where the meson is 
expected to be produced close to where the beam particle 
entered the nucleus. 

The model treats the photon-nucleus reactions as a 
two-step process. In the first step, the incoming photons 
react with a single nucleon taking into account the effects 
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FIG. 8: (color online) Normalized combinatorial background 
(blue) for individual targets compared to data (black). 

of shadowing. In the second step, the produced parti- 
cles are propagated explicitly through the nucleus allow- 
ing for final-state interactions. This step is governed by 
the semi-classical BUU transport equations. A complete 
treatment of the e + e~ pair production from jA reactions 
at Jefferson Lab energies can be found in Ref. (4?J . The 
output is the e + e _ mass spectrum from 7 decays: the 
direct vector meson decays p — > e + e~, lu — > e + e~, and 
<f> — > e + e~, as well as 4 Dalitz decays: A — >■ Ne + e~, 
rj — » e + e~7, 7T° — » e + e _ 7, and ui —> n e + e~ . 

The simulations were employed in two ways. The 
first application was to calculate the detector acceptance. 
Fig. [9] shows the e + e~ detection acceptance for the p de- 
cay as a function of the invariant mass for carbon targets 
at different ^-locations in CLAS. The smooth, slowly- 
varying acceptance does not distort the shape of the in- 
variant mass spectrum. 

The second application was to remove the u>- and <f>- 
meson contributions from the total invariant mass spec- 
trum for each target. The u>- and <p- mesons have long 
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FIG. 9: CLAS acceptance as a function of the e + e~ invariant 
mass for C targets at different z-locations in CLAS. The indi- 
vidual targets are located at z — —12.0, —7.0, 2.0, andS.O cm. 
The results are shown for the p channel and summed over all 
photon energies. 



lifetimes (or = 23.4 fm and 44.4 fm, respectively [39() and 
momenta greater than 0.8 GeV; therefore, their proba- 
bilities of decaying inside the nucleus are low. The to and 
4> channels in the simulation have been treated for many- 
body nuclear effects as they escape the nucleus, while 
changes due to chiral symmetry restoration have not been 
included. A comparison of the simulated lu- and 0-mcson 
line shapes without and with the modifications of Ref. [2] 
shows negligible differences. The lu mass shape contained 
both the direct and the Dalitz decay channels. These two 
channels were combined by fixing the branching ratios in 
the BUU generator. A single w-meson distribution was 
used in the global fit. The acceptance-corrected BUU 
mass shapes for the p, to, and <p mesons were scaled to 
match the experimental mass spectra. Fig. [TU] shows the 
result of the fits to the mass distributions from the 2 H, C, 
and Fe — Ti data after the combinatorial background was 
subtracted. The dot-dashed, dashed, and dotted curves 
are the fits for the p, lu, and <f> mesons, respectively. The 
solid curve is the sum of the three fit contributions. The 
normalized lu- and 0- meson shapes were subtracted. A 
substantial contribution from the p meson is present in 
our mass spectra unlike in the KEK analysis (see Fig. 1 
in Ref. (H). 
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FIG. 10: (color online) Result of the fits to the e + e~ invari- 
ant mass spectrum obtained for the 2 H (top), C (middle), 
and Fe — Ti (bottom) data. The curves are Monte-Carlo cal- 
culations by the BUU model 48, 49] for various vector meson 
decay channels. 



B. The p Mass Spectra 

After subtracting the combinatorial background and 
the lu- and 0- meson contributions, all that remains in 
the e + e _ invariant mass spectra were the p mass distri- 
butions (see Fig. [TTj) for the various targets. The exact 
functional form for the mass spectrum has been obtained 
by calculating the cross section for p-meson production 



including the leptonic decay width [49J and is given by 
2 p 2 T(p) 



7T {p? - m p 2 ) 2 + ^l 2 ^ 2 (^l) , 



(ii) 



where T(fi) is the width of the resonance, M p is the p 
pole mass, and p, is the e + e _ invariant mass. A good ap- 
proximation for the function in Eq. [TT] is a Breit-Wigncr 



10 



function divided by p 3 . In the Vector Dominance Model, 
the photon propagator has the form of 1/q 2 = 1/p 2 - 
This term contributes as l//x 4 in the cross section, while 
the phase space also gives a factor of p, resulting in a 
1/p 3 factor [13, H3, Sll- Indeed the fits to the Breit- 
Wigner//! 3 , rather than a simple Breit-Wigner function, 
describe the data very well (see Fig. [T2|) . For exam- 
ple, for the 2 H target, a simple Breit-Wigner fit gives 
a x 2 P er number of degrees of freedom (Xndf) ec l ua l to 
3.9, while a Breit-Wigner / p 3 gives a Xndf = 1-08 ( see 
Fig. [T2|) . A similar effect is seen with the C and Fe — Ti 
data. For the C spectrum, the fits without and with 
3 f Q „w nnA,,™ „2 of 2.92 and 1.23, 



the l/p 3 factor produce Xndf 



respec- 



tively. For the Fe - Ti data, the Xndf are 302 and !- 35 ' 
respectively. The sensitivity of the fits to the l/p 3 factor 
indicates that the systematic uncertainties in the back- 
ground subtraction procedure are insignificant. Similar 
results are obtained for the heavier targets C and Fe — Ti, 
with larger uncorrelated background (proportional to the 
atomic number of the target). 



VI. RESULTS 

Ratios were made of the mass distributions from the C 
and Fe — Ti data to the 2 H data. Due to the statistical 
level of the data, making the ratios was more sensitive 
to mass and width changes than fits to individual mass 
spectra. Furthermore, taking the ratio of mass spectra 
from different targets minimizes systematic effects. 



A. Ratio of Mass Spectra 

A simple and intuitive way of estimating the sensitiv- 
ity to a mass shift is to measure the slope of the ratios. 
The ratios of the data are shown in the bottom panels 
of Figs. [T3] and [T3J A linear fit was performed to ap- 
proximate the shape of the ratio of the two mass distri- 
butions between 0.55 GeV and 0.80 GeV, and was used 
to estimate the sensitivity to a mass shift. In order to 
illustrate this sensitivity, the ratios were also obtained 
using the simulation with five times more statistics than 
in the data. The result is shown in the top and middle 
panels of Figs. [H (C target) and QJ] (Fe - Ti targets), 
without and with the mass-shift prediction of Hatsuda 
and Lee @, respectively. 

For the extraction of the a parameter, only the analy- 
sis with the Fe — Ti data, where a mass shift is predicted, 
will be discussed. The fit to the data (fT4l bottom) shows 
a slope of 0.46 ± 0.35. In Fig. Qj)J the slope of the lin- 
ear fit to the ratio is shown as a function of the a mass 
shift parameter for a Fe nucleus. The line was obtained 
from simulation by shifting the p mass distributions of 
the Fe — Ti targets and fitting the ratio of the shifted 
to unshiftcd distributions by a polynomial of order one. 
The shaded box on Fig. [15] is our result from the data 
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FIG. 11: Individual Breit-Wigner//^ 3 fits to the p mass spec- 
tra for 2 H (top), C (middle), and Fe - Ti (bottom) 



and gives an a = 0.02 ± 0.02 or |Am| = 10 ± 10 MeV 
(consistent with zero). The prediction of Hatsuda and 
Lee [2j tranlates to a mass shift of 80 ± 30 MeV 

Our results, along with the prediction of Hatsuda and 
Lee |2j and the result of KEK [30( , are shown in Fig. [THl 
A comparison to the theory and other published results 
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FIG. 12: Fits to the p mass spectrum for the deuterium 
data using a Breit-Wigner (dashed) and a Breit-Wigner//i 3 
function (solid). 
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FIG. 14: (color online) Ratios of the e + e~ invariant mass 
for Fe — Ti to 2 H. An explanation of the fits is given in the 
caption of Fig. 1131 



is given in Sec. IVIII1 

The result of the fits to the ratio obtained from the 
simulations with no shifts (Figs. [T3l and [Ml top), indi- 
cated essentially the same result when compared to the 
data (Figs.[THland[T4l bottom). Note that the deuterium 
target in the BUU model was not treated as a nucleus, 
but as a proton and a neutron. Therefore, some nuclear 
effects (mentioned in Sec. IV A[) led to a non-zero slope of 
the ratios. 



B. Simultaneous Fits to Mass Spectra and Ratios 



FIG. 13: (color online) Ratios of the e + e~ invariant mass for 
C to 2 H. The ratios were obtained from the data (bottom), 
and the simulations with (middle) and without (top) the mass 
shift predicted by Ref. |2§. The number of simulated events 
is five times more than the current data. The slope of a 
simple linear fit shown on each plot is used to investigate 
the sensitivity of the current data to a mass shift (shown in 
Fig. El). 



A more precise method of measuring the mass and 
width of the p meson in various targets was developed. 
These values were obtained by performing a simultaneous 
fit to the mass spectra and the ratio of each spectrum to 
the 2 H data. The simultaneous fits to the mass spectra 
and ratio were performed to impose more constraints on 
the fits. The results of the fits are shown in Figs. [T71 and 
1181 and are summarized in Table |TJ 
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FIG. 15: (color online) Slope (1/GeV) of the linear fit to 
the ratio of the p mass distributions in the Fe — Ti data as a 
function of the mass shift parameter a. The prediction of the 
Hatsuda and Lee model [j| (blue) and the KEK result [3(| 
(red) are compared to the data within ±1<t (shaded green 
box). 
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FIG. 16: The shift parameter a, with the errors given hori- 
zontally, from the slope analysis with the Fe — Ti and 2 H data 
(circle). Also shown is a comparison with the KEK result [3(J ] 
(square) and the prediction of Hatsuda and Lee model [2] 
(triangle). For the KEK result, the error is smaller than the 
symbol. 



VII. SYSTEMATIC UNCERTAINTIES 

The major sources of systematic uncertainties are sum- 
marized in Table [ill To estimate the systematic uncer- 
tainties, the measured or nominal value of the shift pa- 
rameter a was compared to the a values determined when 
using alternative cuts or corrections. The a parameter 
was determined by the fit to the linear slope of the ra- 
tio of mass distributions as described in Sec. IVI Al The 
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FIG. 17: (color online) The result of a simultaneous fit to 
the p mass spectra for 2 H (top), C (middle), and the ratio of 
C/ 2 H (bottom) 



difference in a was utilized as a measure of the system- 
atic uncertainty. The estimated common uncertainty for 
a given effect was the weighted root-mean-square of the 
difference in a for all the altered cuts or corrections. For 
every case except the target position study, the Fe — Ti 
and 2 H data were analyzed. 

For each vertex cut listed in Table HH loose and tight 
cuts around the nominal value were applied. 

For the momentum corrections and target energy loss 
corrections, the values of a were compared with and with- 
out the corrections. The target position uncertainty was 
obtained by comparing the invariant mass spectrum from 
the first two carbon targets (z = —12.0 cm and —7.0 cm) 
to the last two carbon targets (z = —2.0 cm and 3.0 cm). 
The Aa was about 0.006 due to the low statistics of the 
same-charge events that contributed to the estimation of 
combinatorial background for each carbon target. The 
total systematic uncertainty is much smaller than the 
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FIG. 18: (color online) The result of a simultaneous fit to the 
p mass spectra for 2 H (top), Fe — Ti (middle), and the ratio 
of Fe-Ti/ 2 H (bottom) 



statistical uncertainty. 

The systematic uncertainty in the background sub- 
traction was evaluated with two fits to the experimen- 
tal p mass spectra, a Brcit-Wigner shape and a Brcit- 
Wigner//i 3 . As stated in Sec. IV Bl the systematic un- 
certainty due to the shape of the background function is 
estimated to be negligible. For the background normal- 
ization factor, the 7% statistical uncertainty from the 
mixed-event technique (see Sec. IIVC|) was propagated 
through the simple ratio analysis to estimate Act. The 
result is a 0.001 difference in the measured a. 

As discussed before, the structure seen in the ratio 
plots can be very sensitive to a mass shift and/or change 
in the width. More realistic functional forms describing 
these changes were studied, but a linear fit was used as 
a simple measure for observing small possible changes 
in the mass. The characteristic plot of the ratio of two 
spectra indicated the sensitive range to the mass shift. 



TABLE I: The mass and width of the p meson obtained from 
the simultaneous fits (M exp ,T exp ) to the mass spectra for each 
target and the ratio to 2 H compared to the result of the 
simulations (M 3im ,T sim ). The masses and widths are con- 
sistent with the natural values [39| ] (770.0 ± 0.8 MeV and 
150.7 ± 1.1 MeV, respectively) adjusted for the collisional 
broadening. The units are in MeV. 



Target 






M ■ 

lvl sim 


r 


2 H 


770.3 ±3.2 


185.2 ±8.6 


774.5 ± 4.9 


160.1 ± 10.2 


C 


762.5 ±3.7 


176.4 ±9.5 


773.8 ± 0.9 


177.6 ±2.1 


Fe,Ti 


779.0 ±5.7 


217.7 ± 14.5 


773.8 ± 5.4 


202.3 ± 11.6 



TABLE II: Summary of systematic uncertainties in the mass 
shift parameter a. For the total error, the individual contri- 
butions were added in quadrature. 



Source 


Description 


Aa 


Vertex Cuts 


Radial Position 


0.003 




z Position 


0.001 




Time 


0.002 


Corrections 


Momentum Corrections 


0.0015 




Target Energy Loss 


0.001 


Target Position 


z Dependence 


0.006 


Background Sub. 


Mixed- Event 
Normalization Factor 


0.001 


Fit Range 


Table [nil 


0.007 


Total 


Sum of individual sources 


0.01 



Using simulation data for known mass shifts, the slope 
of the linear fit was translated to a and compared with 
the slope in Fig [151 The systematic uncertainty associ- 
ated with the choice of the fit range is summarized in the 
Table Mil For the determination of the systematic un- 
certainty, the fit range was varied by 0.02 GeV with all 
of the permutations taken into account (varying only the 
upper limit , varying only the lower limit, and varying 
both limits). The total value is an over-estimation of the 
systematic uncertainty due to the choice of the fit range. 



TABLE III: Summary of the results of the systematic uncer- 
tainty on a due to the choice of fit range. The change to the 
limits in each case was 0.02 GeV. 



Fit Range 


Sys. Uncertainty in a 


Increase/decrease upper limit 


0.003 


Increase/decrease lower limit 


< 0.001 


Increase/decrease both limits 


0.006 


Shift the fit range 


0.003 


Total 


0.007 
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VIII. DISCUSSION AND CONCLUSIONS 

This experiment successfully detected the light vector 
mesons via their rare decay into e + e~ pairs, in order 
to eliminate final-state interactions. The CLAS detec- 
tor is ideal for this measurement as it can discriminate 
effectively between lepton and pion pairs to the level of 
10~ 7 . Background contributions from Bethe-Heitler pro- 
duction, 2ir° Dalitz decays, and combinatorial processes 
were removed. The determination of the combinatorial 
background was possible due to e + e + and e~e~ samples 
in the data. The narrow lu- and 0-meson peaks were re- 
moved. What remained in the mass spectra was the ex- 
perimental p-meson distribution described very well by 
a Breit- Wigner distribution scaled by 1 //i 3 (see Eq. [TT]) . 
This experiment had the unique characteristics of an elec- 
tromagnetic interaction in both the production and de- 
cay of the vector mesons. The p-meson mass spectra 
have been extracted for 2 H, C, and Fe — Ti nuclei in a 
model-independent way. With their long lifetimes and 
momenta greater than 0.8 GeV, most lu and ip mesons 
decay outside the nucleus and were treated as in- vacuum 
decays. 

By analyzing the ratio of the Fe — Ti to the 2 H mass 
distributions, a value of 0.02 ± 0.02 for the mass shift 
parameter a was obtained. 

This result differs from the KEK measurement [30| . 
The measurement at KEK detected the p, lu, and <p 
mesons from the e + e~ decay. Unlike the present exper- 
iment, the vector mesons were produced with a proton 
beam at an energy of 12 GeV. The targets in the KEK 
experiment were C and Cu. The C data was taken as 
the reference where no medium modifications were ex- 
pected due to its small effective density. The lu- and 4>- 
mesons were visible in the e + e~ invariant mass spectra. 
The KEK analysis did not have a sample of same-charged 
leptons by which to extract the normalization of the com- 
binatorial background. Instead, the background contri- 
bution was fit along with the lu- and </>-meson shapes. 
Without an absolute determination of the combinatorial 
background, the p-meson signal was suppressed and in- 
cluded in the background shape. The lost p-meson yield 
was recovered by using a theoretical model to constrain 
the p meson to lu meson production strengths. Their 
background subtraction procedure lead to the question- 
able result of a = 0.092 ± 0.002. 

In a comparison with theoretical predictions, the orig- 
inal calculations of mass shifts by Refs. Q and are 
ruled out by our result under our experimental condi- 
tions. Predictions by the Giessen group [2(| and the Va- 
lencia group [1H are consistent with our result of a small 
to no mass shift. 

Beyond the mass shifts, there can be changes to the p- 
meson width. From our analysis of the simultaneous fits 
to the mass spectra and their ratios, the extracted widths 
are consistent with collisional broadening and show no 
signs of further modifications. In the KEK analysis, the 
width of the p-meson was fixed. Again, our result is 



consistent with Refs. [20l [2l|. 

Besides the investigations of the p-meson properties, 
the yields of the lo and <f> mesons can be extracted. From 
Fig. [TU1 the strength of the lo and (j> signals decrease as 
the target mass increases as a result of absorption. An 
analysis of the nuclear transparencies is underway to ac- 
cess the in-medium widths of these mesons. Substantial 
increases to the in-medium widths of these mesons have 
been predicted by Refs. [H, M, SI US US 113, M, El . 

From the p mass spectrum for p p > 0.8 GeV, there is 
no evidence of large modifications to the mass and width. 
Our analysis with the a parameterization is compatible 
with no mass shift. We set an upper limit of a = 0.053 
with a 95% confidence level. 

The next step is to extract the in-medium spectral 
function (60| . However, obtaining the spectral function 
from the mass spectrum is not trivial. Pointed out in 
Ref. (6lj |. the mass spectrum contains information about 
the spectral functions, partial decay widths, and the 
coupling constants, all of which can be modified in the 
medium. In light of our result, more theoretical work is 
needed in extracting the spectral functions. 
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